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Abstract 

Short title: Non-classical plasticity theories and generalized normality. 

A simple way to define the flow rules of plasticity models is the assumption of 
generalized normality associated with a suitable pseudo-potential function. This 
approach, however, is not usually employed to formulate endochronic theory and 
non-linear kinematic (NLK) hardening rules as well as generalized plasticity mod- 
els. In this paper, generalized normality is used to give a new formulation of these 
classes of models. As a result, a suited pseudo-potential is introduced for endochronic 
models and a non-standard description of NLK hardening and generalized plasticity 
models is also provided. This new formulation allows for an effective investigation 
of the relationships between these three classes of plasticity models. 
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1 Introduction 



The models proposed so far in the literature to describe the rate independent 
inelastic behavior of real materials subjected to monotonic or cyclic loading 
conditions can be essentially classified into two main families: (i) models where 
the present state depends on the present value of observable variables (total 
strain, temperature) and of suitable internal variables; (ii) models, indicated 
here as hereditary, that require the knowledge of the whole past history of 
observable variables. 



The first group encompasse s, for instance, the classi c al mo dels of Prandtl- 
Reuss and Prager (see e.g. iLemaitre and Chabochd (119901 )) and the NLK 
hardening model of lArmstrong and Frederick! (119661 ). in its origi n al for m as 
well as in the m o dified versions recently proposed by IChabochd (119911 ) and 



Ohno and Wand (119931 ) in order to improve the ratchetting modelling. For 
these models, the well known notions of elastic domain and loading (or yield- 
ing) surface apply. Associativity and non-associativity of the plastic strain 
flow rule are also well-established concepts, as well as the assumption of 
generalized associativity (or generalized normali ty), relating all internal vari- 
able flow directions to a g iven loading surface (IHalphen and Nguvenl. I1975T) 
( Jirasek and Bazantl . l2002l ). Using the language of convex analysis (iRockafellarl . 



19691 ). generalized normality entails that the flows of all internal variables be- 



long to the subdi fferentia l set o f a given s c alar n on-negative function called 
pseudo-potential ( jMoreaul . Il970l ) (IFremondl . 120021 ). 



Among internal variable theories, generalized plasticity deserves special at- 
te ntion. A first important step for its formulation was the idea, suggested 
by lEisenberg and Phillips! (Il97ll ). of a plasticity model where, despite classi- 
cal plasticity, loading and yielding surfaces are not coincident. Then, start- 
ing from an axiomatic approach to describe inelastic behavior of materials, 
Lubliner proposed some simple generalized plasticity models, able to represent 
some observed experi mental behavior of metals (Lubliner, 1974, 1980, 1984) 
( jLubliner et all Il9931 ). More recently, general ized plasticity has been used f or 
describing the shape memory alloy behavior ( Lubliner and Auricchio . 19961 ). 



Endochroni c models (jValania . Il97ll ) and Bouc-Wen type models (lBoud . ll97ll ) 
( jWenl . Il976l ) are two important examples of hereditary models. Endochronic 
theory has been developed during the seven ties and used f or modelling the 

plastic behavior of metals (see, for instance, Valanid ( 1971 ). Valanis and Wu 

( 19751) ') and the inelastic b e havio r of concrete and soils (among others, Bazant and Krizek 
(ll976h . iBazant and Bathl (Il976h . The endochronic stress evolution rule de- 



pends on the so-called intrinsic time and is formulated by a convolution inte- 
gral between the strain tensor and a scalar function of the intrinsic time called 
memory kernel. When the kernel is an exponential function, an incremental 
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form of endochronic flow rules exists, which is commonly used in standard 
analyses and applications. 



Models of Bouc-Wen type are widely e mployed for modelling t he cyclic behav- 

ior of structures in seismic engineering (IBaber and Wenl . ll98ll ) (jSivaselvan and Reinhornl . 



200CH ) and for representing h ysteresis of magn et o-rheological dampers in semi- 



active control applications (jSain et al.l . 119971 ) (IJansen and Dykd . 120001 ) . The 



strict relationship between endochronic and Bouc-Wen typ e models has been 
menti oned several times in the literatu re (see, among other s . iKarrav and Bouc 
fll989h and lCasciatil fjl989h ). Recently. lErlicher and Pomtl (l2004 ) showed that 
the fundamental element of this relationship is the choice of an appropriate 
intrinsic time. 



Endochronic theory and classi cal inter n al var iable theory have been compared 
by using several approaches: iBazantl (119781 ) observed that for endochronic 
theory the notion of loading surfa ce ca n still be introduced, but i t looses 
its physical meaning; IValanisI (119801 ) and IWatanabe and Atluril (119861 ) proved 
that a NLK hardening model can be derived from an endochronic model by 
imposing a special intrinsic time definition. Moreover, a comparative study 
be tween NLK hardening and g eneralized plasticity models has been presented 
by lAuricchio and Taylorl (119951 ). A tight relationship between endochronic the- 
ory and generalized plasticity is also expected to exist, but, by the authors' 
knowledge, no analysis on this subject exists. More generally, there is a lack 
of unified theoretical framework, on which formal comparisons between these 
plasticity theories could be based. The main goal of this paper is the formu- 
lation of t his theoretical fr amework using the classical notion of generalized 
normality (Moreau, 1970h ( Halphen and Nguyen . 1975 ). As a result, a new 
formulation of endochronic and NLK hardening models as well as general- 
ized plasticity models is suggested and is used to investigate the relationships 
between them. 



The paper is organized as follows: in the first section, the standard theoret- 
ical framework of thermomechanics is briefly recalled, with reference to the 
notions of pseudo-potential and generalized normality as well as Legendre- 
Fenchel transform and dual pseudo-potential. In the following sections, several 
plasticity models are presented and are shown to fulfil the generalized normal- 
ity assumption. The Prandtl-Reuss and endochronic models are considered 
first, in both standard and multi-layer formulations. Then, NLK hardening 
model and generalized plasticity follow. The discussion is limited to initially 
isotropic materials, whose plastic behavior is governed by the second invariant 
of the deviatoric stress, J2, known as von Mises or J2 materials. No stability 
analysis is provided, as it is beyond the purposes of this contribution. 
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2 General thermodynamic framework 



Under the assumption of infinitesimal transformations, the classical expression 
of the local form of the first and second principle of thermodynamics can be 
written as follows: 



Ts 



— sT — div (q) + r + cr : i 



<f> s (t) = s + div - 1 > 



(1) 
(2) 



where the superposed dot indicates the time derivative; s is the entropy density 
per unit volume, q is the vector of the flowing out heat flux, T is the absolute 
temperature, r is the rate of heat received by the unit volume of the system 
from the exterior; cr is the second order symmetric Cauchy stress tensor; e 
is the tensor of small total strains; $ s (t) is the rate of interior entropy pro- 
duction. In the vector space of all second order tensors, the Euclidean scalar 
product : is defined by x : y = x^y^ the vector subspace of second order 
symmetric tensors is denoted by S> 2 . The Helmholtz free energy density per 
unit volume is a state function defined as 

* = *(e,T, Xl , = *(v) (3) 



where Xn ■■■iXni are the tensorial and/or scalar internal variables, related to 
the non-elastic evolution and v = {e, T, Xn •••> Xn} * s the vector containing all 
the state variables, namely the total strain tensor, the temperature and the 
internal variables. 

For isothermal conditions, the use of Eq. ( [[]) in the inequality (j5J) leads to 
$ m (t) := T$ s (t) = Ts = o- : i - * > (4) 



which states that the intrinsic or mechanical dissipation $ m (rate of energy 
per unit volume) must be non-negative. The non-dissipative thermodynamic 
forces are define d as functions of the free energy density ^ (see, among others, 
Fremondl (12002^ 



ds ' 



-nd ._ 



dXi 



q 



nd 



9* 
dv 



(5) 



The non-dissipative stress cr is associated with the observable variable e, 
while r" d are associated with the internal variables x%- All non-dissipative 
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forces can be collected in q ad = [cr nd , T™ d , r^f j. Hence, by substituting ((SJ) 
into (prj, one obtains 

JV 

$ m (t) = (<r - : e - 2 < d ■ X l = * ■■ e - q nd ■ v > (6) 

i=i 



where v is the vector of the fluxes, belonging to a suitable vector space V. The 
vector spaces considered in this paper are isomorph to M. n and the same hold 
for their dual V* (see Appendix). The symbol ■ indicates the scalar product 
of two objects having the same structure: two tensors, two scalar variables 
or two collections of tensorial and/or s c alar y ariables (the same notation has 



been used e.g. by iJirasek and Bazantl (120021 . pg. 428)). The inequality 



can be written in a slightly different manner, by introducing the dissipative 
thermodynamic forces 



.d 

1) 



T N 



nd 



_nd 
" T l ) 



rid 



e V* 



(7) 



Hence, 



N 



i=i 



The forces q d have to be defined in such a way that the couple (q d , vj always 
fulfils the inequality (jHJ). Therefore, some additional complementary rules have 
to be introduced. They can be defined by assuming the existence of a non- 
negative, proper, convex and lower semi-continuous function : V — ► (— oo, oo] 
(Appendix, item 2), called pseudo-potential, in general non-differentiable, such 
that <p (0) = and: 

q d G^(v) (9) 



where d indicates the sub-differential operator (see Appendix, item 3). This 
condition is called generalized normality. A more detailed way of writing (Q 
is 

q d e <90(v';v)|., = . (10) 



As a matter of fact, is a general function of the fluxes v' and may also 
depend on the state variables v. However, the sub differential is taken, by 
definition, only with respect to the fluxes v' and the thermodynamic force q d 
corresponds to the subdifferential of at v' = v, where v is the actual flow. By 
using the properties of sub-differentials, it can be proved that for dissipative 
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forces denned by (jHJ), the inequality q d ■ v > is always fulfilled (Appendix, 
item 4). Therefore, the second principle (jSj) is also satisfied. 

A dual pseudo-potential 0* : V* — > (—00, 00] can be defined by the Legendre- 
Fenchel transform of 6: 



0* (q d ') 



sup 

v'ev 



q d ' • v' - (v')) (11) 



When has an additional dependence on the state variables v , then ffTTT) 
leads to 0* = 0* (q d '; v). It can be proved that the dual pseudo-potential is a 

non- negative, proper, convex and lower semi- continuous function of q d ' , such 
that <fi* (0) = (see Appendix, item 5). The dual normality condition reads 



v G 



d<P* (q d ) (12) 



where q d is the actual value of the dissipative force. The expression ( TT21 is 
equivalent to 



v G 



(<i^)Lw < 13 > 



and it guarantees that q d • v > (Appendix, item 5). M oreover, it defines the 



complementarity rules of generalized standard materials (Halp hen and Nguyen 



1975|), sometimes called fully associated materials (jJirasek and Bazantl . 12002 



pg. 452). 



Plasticity is characterized by a rate-independent memory effect (iVisintinl . 119941 . 
pg. 13). This special behavior occurs when the pseudo-potential is a posi- 
tively homogeneous function of order 1 with respect to the fluxes v'. In this 
case, provided that q d is computed from (J9j) or that v derives from (fl2|) . it can 
be proved that the pseudo-potential at v is equal to the intrinsic dissipation, 
viz. (v) = q d ■ v = $ m (Appendix, item 6). Moreover, the dual pseudo- 
potential 0* becomes the indicator function of a closed convex set EcV* and 
the normality rule ffl2l entails that, given the dissipative force q d G E, the 
flux v fulfils the following condition: 

Vq d ' G E (q d> - q d ) ■ v < (14) 



viz. for a given dissipative force q d , the flow v defined by ( 1141) (or, equiv- 
alently, by (|T2l) or (fTBl ) is such that its power when it is associated to the 
actual force q d is always greater or equal to the power q d • v of all the 
othe r dissipative forces q d G E (generalized maximum-dissipation princi- 
ple ( Halphen and Nguyenl . 19751 )). When q d G <9E, the inequality (TT41) in- 
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dicates that v belongs to the cone orthogonal to <9E at the point q d . When 
q d G int (El, it forces the flow v to be zero (Appendix, item 7). 



3 Prandtl-Reuss model 

3. 1 Perfectly plastic Prandtl-Reuss model 

In order to illustrate the general procedure that is adopted hereinafter, the 
basic example of the Prandtl-Reuss model is considered first. The relevant 
state variables are the total and the plastic strain v = (e, e p ) and q nd = 
(cr nd , T nd ^j are the associated non-dissipative thermodynamic forces. The usual 
quadratic form of the Helmholtz free energy density \I/ is used, in order to 
preserve the linear dependence of all non-dissipative forces with respect to 
state variables: 

* = - (e - e p ) : C : (e - e p ) (15) 

For isotropic materials, the fourth-order tensor of the elastic moduli is equal 
to C = (K — |G) 1 (g> 1+2GI, where K is the (isothermal) bulk modulus, G is 
the shear modulus and ® is the direct (or outer) product of two second order 
tensors. The assumption of isotropy is always adopted, even if the concise 
symbol C is used. The non-dissipative forces associated with (fl5l) can be 
derived by means of 

<r nd = C : (e - e p ) , r nd = -C : (e - e p ) (16) 

The evolution of the dissipative forces q d = (er d , r d ^ G S 2 x S 2 := V* is defined 
by introducing a suitable pseudo-potential 0, which is a function of the fluxes 
v' = (i 1 , i p "j e S 2 x S 2 := V (x is the cartesian product): 

D = i p ') G V such that tr (i p ') = o} 

where tr(u) indicates the trace of the tensor u G§ 2 . The norm of the second 
order symmetric tensor u is given by ||u|| = y/Uij Uij. If in addition tr (u) = 
Ua = 0, then ||u|| 2 = 2 J 2 (u) where J 2 (u) is the second invariant of the 
deviatoric part of u; a y is the one-dimensional tension stress limit and % 
is the indicator function of the set D, namely % = if tr (i p J = and 



e p 
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Ig = +00 elsewhere. This set is the effective domain of (Appendix, item 
2). The pseudo-potential is a homogeneous function of order 1 with respect 
to (i 1 , i p J and therefore a rate independent constitutive behavior is expected 

and the dissipation $ m is equal to ^/|o" 3/ ||e p ||, where i p is the actual plastic 
flow (Appendix, item 6). The indicator function % accounts for the fact that 
plastic deformation occurs without volume changes (plastic incompressibility) . 
This assumption is usual for metals and has been validated by experimental 
evidence. 

The pseudo-potential (ft*, dual of 0, can be computed using the Legendre- 
Fenchel transform (Appendix, item 5) and is equal to: 

0* (>', r*') = su P( ,, , P>S (** : a + : e' - 0) 



The dual pseudo-potential 0* is the indicator function of a closed convex set 
E. Hence, = 0** is the support function of the same set (Appendix, item 6). 
Moreover, since does not depend on i', the dual pseudo-potential 0* can be 
written as the sum of two indicator functions (Appendix, item 7): 



(a d \r d ')=l Q (>')+I E (r d ') 
G § 2 such that / 



E 



dev (r d ')| - J\ay < 0} 



(19) 



where dev(u) is the deviatoric part of u G§ 2 . The first term is equivalent to 
the condition cr d = 0, while the other indicator function 1^ defines a region 
in the T d stress space. Recalling that the actual value of T d , viz. r d , fulfils 
the condition T d = —r nd and that the only possible value for er d is zero, 
using (fl6l) it is straightforward to see that T d = a = cr nd and that E can 
also be interpreted as a set in the cr stress space. The associated function / is 
known as loading function and the condition / = defines the plastic states. 
The interior of E is associated with the elastic states and the whole (closed) 
set E contains all plastically admissible states. The actual flows (e, i p ) can 
now be derived from (Tl9]) by computing the sub differential set of 0* and then 
considering it at (cr d ',T d '^ = (^cr d J T d ^j. Hence, no restrictive conditions are 
imposed on e, while the plastic strain flow reads (Appendix, item 7) 

e p = JP^Xx = n A with 

IIM T )II (20) 

A > 0, / (r d ) < 0, Xf (r d ) = 
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Observe that / in the loading-unloading conditions in the second row of (120]) is 
computed at the actual stress state. The plastic multiplier A is then evaluated 
by im posing the consistency condition, i.e. Xf = (see e.g. ISimo and Hughes 
jl988h ). with 



df 
dT d> 



:r d ' 



(21) 



Note that / is computed from the general expression of / \T d j and then evalu- 
ated at the actual state T d> = T d . Consistency corresponds to the requirement 
that in order to have A > 0, the actual dissipative force T d G <9E cannot 
leave <9E during the plastic flow. Hence, by using ( 1201) and the relationship 
r d = <t = C : (e — s p ), one has 

dev I T d ) 

/ = ) J : dev (t d ) = n : C : e-n : C : n A (22) 
\\dev(r d )\\ V J K J 



thus X = H{f) = if (/) (n : e) , where (x) = ^ (McCauley brackets) 

and if (/) is the Heaviside function, equal to zero for / < and equal to 1 
elsewhere. 

In summary, the Prandtl-Reuss perfectly plastic model was formulated by 
means of the Helmholtz free energy \l/ and the pseudo-potential <p\ then, the 
dual potential 0* was computed from the Legendre-Fenchel transform of 0; 
the sub differential set of <fi* was used to define the fluxes and the consistency 
assumption led to the determination of the plastic multiplier A. In the next 
sections, this approach will be used to formulate two Prandtl-Reuss models 
with isotropic hardening and other more complex plasticity models, such as 
endochronic, NLK hardening and generalized plasticity models. In order to 
get this result, some non-standard expressions for the pseudo-potentials (j> and 
0* are introduced. 



3.2 Classical Prandtl-Reuss model with isotropic hardening 



The vector of the representative state variables for a Prandtl-Reuss model 
with isotropic hardening is equal to v = (e,e p ,()i while q nd = (^er nd , er , R nd ^j 
are the associated non-dissipative forces. The scalar internal variable ( and the 
associated forces R d and R nd are introduced in order to represent the isotropic 
hardening. Moreover, v = (i, i p ,(j G V = S 2 x S 2 x R is the flux vector and 

q d = (cr d , cr d } R d j G V*= S 2 x S 2 x R contains all dissipative thermodynamic 
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forces. The Helmholtz free energy is assumed to be of the form 



¥ = i(e- e '):C:( e -e*)+e(C) (23) 



where £ (() is a scalar function such that £ (0) = and ^ (0) = 0. It follows 
that 

<r nd = C : (e - e p ) , r nd = -C : (e - s p ) , R nd = % (C) ( 24 ) 



The pseudo-potential is assumed of the following form: 

4>(i',i p ',c) = vik c'+Id (^^',c') 

(25) 

D= | (e', £ p ', C') G V such that tr (>') = and C' > ||e p '|| } 



The first term in the expression of <fi is the same as in the perfectly-plastic 
model when (' is equal to the norm of i p . The second term, that is the 
indicator function %, depends not only on tr(i p ), but also on the flow 
which is forced be greater or equal than the norm of the plastic strain flow. 
This inequality guarantees that (' and <ft are non-negative and entails that D 
is convex and closed (see Appendix, item 1 and Figure [T^i, which illustrates 
the projection of D on the (^i p -plane for the tension-compression case). 
The dual pseudo-potential is different from ([TBI) , due to the presence of the 
dissipative force R d associated with 



(cr d \ r d \ R d ') = sup ( ^ c , )eB (a d ' : i> + r d ' : e p ' + R d ' (> 



(26) 



where E = {(r d ',R d ') G § 2 xE such that / (r d ',R d ') < o} and 




The loading function / defines a convex and closed region E in the \r d ' , R d J 

space, where the actual value of R d> , viz. R d = —R nd = —^p- governs isotropic 
hardening (or softening). The limit stress becomes greater than its initial value 
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\j\(J y when > and less when < 0. Figure [Tb illustrates the set E. 
The flow rules follow from the generalized normality conditions: 



e' J 



devl T d 



-A = nA, 



C = A 



\\dev[T d ) 

with A > 0, / < 0, A/ = 



(28) 



The flow of the internal variable ( is equal to the plastic multiplier A, which 
can be evaluated by imposing the consistency condition: 



/ 



% : f + SL R« 



dr 



dR d ' 







. (T d '=T d ,R d '=R d ) 



(29) 



It follows that 



(n : C : e) 



n : C : n-F 



H(f) 



(n:i) 



dC 2 



1+ 



1 aPg(C) 



(30) 



2G d( 2 



where H (f) and ( ) still indicate the Heaviside function and McCauley brack- 
ets, respectively. 



3.3 Modified Prandtl-Reuss model with isotropic hardening 



The classical model of the previous section can be extended as follows. Assume 
the state variables v = (e, e p , () and let the Helmholtz energy be equal to 



* = -{e-e*):C:(e-e?)+t{C) 



(31) 



As a result, the non-dissipative forces are the same as in Eq. 
generalized definition of the pseudo-potential is adopted: 



Then, a 



(e\ C; C) = Ulvyg (0 - ^) C' + % i p ', C 

= I (i', i p ', C) e V such that tr (>') = and £' > \\i P '\\ } 



(32) 



In this case, explicitly depends on the internal variable £, by means of 



and of the function g((), positive and such that g (0) = 1. In the particular 
case where g (£) = 1 + ^p-y/f^K the classical expression given in Eq. fl25l) is 
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recovered. The dual pseudo-potential 0* can be evaluated from the standard 
procedure, thus yielding: 



(VV,i? d ';C)=Io (<r d ')+I E (T d \R d ';() 



(33) 



where E = {(r d ', R d ') G § 2 xR such that / (r d ',R d ';() < o} and 



f(T d \R d ';C) = \\dev(T d ')\\-^a y g(C) 



R d ' 



(34) 



In Figured the projection of D on the (s p , -plane and the set E are depicted 
for the tension-compression case, with the assumption £ (() = 0. The flow 
rules are the same as in the previous case and they are reported below for 
completeness: 



.„ dev(T d ) ■ ■ 

£ = n , )_A n A = nA 



C = A 



<feu(T d J 

with A > 0, / < 0, Xf = 



(35) 



In this case, / has to be computed accounting for the state variables. Hence 
consistency condition reads 



/ 



1L . +* + 1L & + ?Lt 

dr d ' ' + dR d ' <9C 



(36) 



=T d .Ft d '=R d ) 



and the plastic multiplier becomes equal to: 

n ■ ^ ■ n +V 3^9 dC i + V3 2G d( 



provided that 1 + y/ff^^ 1 
which occurs when ^jp- < 0. 



> 0. This condition does not prevent softening, 



The comparison of Eqs. (127)1 and (134)) proves to be very interesting. First, the 
usual loading function only depends on the dissipative forces, while / in Eq. 
is also related to the internal variable (. Moreover, since R d = —R nd = 
'^dp~> ^he loading function (1341) at (r d ,R d ^J becomes 



/ (t* R d ; C) 



dev [t 



~ \ n°y9 (0 



(3f 
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This expression shows that the actual limit stress is equal to ^J^a y g(() and 
is independent from the function £ (() introduced in the Helmholtz energy 
density (this is not the case for the classical Prandtl-Reuss model). 

The difference between the two Prandtl-Reuss models can be also explained 
in terms of mechanical dissipation $ m . For the modified Prandtl-Reuss model, 
it is equal to 

*^(/|^(0-^)c (39) 



which is non-negative provided that < y§cr y g (C)- The case of a mono- 
dimensional monotonic loading is depicted in Figure [3) The standard Prandtl- 
Reuss model is characterized by the fact that the energy R d ( associated to 
isotropic hardening is not dissi pated. For this reason it is some times referred as 
energy blocked in dislocations (ILemaitre and Chabochel . ll990l . pg. 402). Hence, 
the mechanical dissipation is equal to yfcr^C for any function £(C)- Conversely, 
for the modified Prandtl-Reuss model the amount of mechanical dissipation 
depends, for a given function g(C), on the choice of £ (£). Figure [3b reports 
the case of generic functions g(() and £(C)- Figure [3b and EH correspond to 
£(C) = and to the case where the modified model is equal to the classical 
one, respectively. 



3.4 Multi-layer models of Prandtl-Reuss type 



Modified Prandtl-Reuss models , defined by Eq s. fl3"TT) - fl3"2l . can be directly 
extended to multi- layer models (jBesselingi . Il958l ). They consist of a system of 
N elastoplastic elements connected in parallel. When every individual elements 
are Prandtl-Reuss models, the corresponding multi-layer model is indicated as 
of the Prandtl-Reuss type. This is the case in the present section. Hence, let 



N 



8=1 



N 

E 

i=l 



■):C: 



(40) 



be the Helmholtz energy density, defined as the sum of N expressions of the 
type f[3"T|) . The internal variable e\ is the plastic strain of the generic element 
i, while Q is the scalar variable associated with the isotropic hardening of the 
same element. All elements have by definition the same elastic modulus ten- 
sor, chosen to be equal to C = jj (K — 1 <8> 1+2GT . The non-dissipative 
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thermodynamic forces read: 



-C:(e-eJ), R t 



nd _ dtjiXCi) 



Let us introduce the pseudo-potential as the sum of iV independent functions 
of the type (1521) : 



{ (e',ef,Ci) E V such that tr (if) = and Q > \\if\\ } 



(42) 



The limit stresses a yi as well as the isotropic hardening functions 9i (Ci) are, 
in general, distinct. The conjugated pseudo-potential is in turn the sum of iV 
independent functions, i.e. 4>* = Y^iLi 4>i with 



if (cr d \ rf, Rf) = sup, , x : e' + rf : if + Rf' Q 

= IoK)+I El (rf,< 



where Ej = {(rf , i?f ) G § 2 xM such that /< (rf , i^; Ci) < o} and 



i' , <*&(Cf 



/* (rf , i?f ; Ci) = \\dev (rf) \\ - ^a yi Qi (Ci) + JSf + 



(44) 



Therefore, iV independent loading surfaces have been defined. Using the stan- 
dard procedure based on the normality assumption, N pairs of flow rules of 
the type (1351) can be derived: 



•p _ dev { T t) X _ { 
Si - \\dev(Tf)\\ Ai ~ UiAi > 



Ci — \ 



with \i>0, fi<0, \ifi = 



(45) 



Moreover, by imposing the consistency conditions and accounting for Eqs. (1411) 
as well as the identities rf = — r" d and Rf = —R^ d , each plastic multiplier 
can be easily determined by an expression of the type ([3] 



Xi = H(f i [ 



(n, : C : e) 



n, : : C : n, 



2 dgjCCj) 
3°V dQ 



H{fi. 



1 -|_ /2 rfgife) 
3 2G dCi 



(46) 
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provided that 1 + 



2 o-yi dgtiCi) ^ n 

3 2G dd U- 



The Distributed Element Model (jlwanl . Il966l ) (jChiang and Beckl . I1994J ) is re- 
covered when gi (Q) = 1 and & (Q) = 0. 



4 Endochronic theory 



Endo chronic theory was first formulated by I Valanid (119711 ). who suggested the 
use of a positive scalar variable called intrinsic time, in the definition of 
constitutive laws of plasticity models. The evolution laws are described by 
convolution integrals involving past values of the state variable e and suitable 
scalar functions depending on d called memory kernels. When the memory 
kernel is exponential, the integral expressions can be rewritten as simple differ- 
ential equations, which, for an initially isotropic endochronic material fulfilling 
the plastic incompressibility assumption, read: 



tr (<x) = 3K tr (e) 

dev (&) = 2G dev (i) —(3 dev (tr) $ 



(47) 



with (3 > 0. These relationships are equivalent to: 
cr= C : (e - e p ) , 

< C = (K - f G) 1 ® 1+2GI , (48) 
tr (e p ) = and i p = $ 



where •& > is the time-derivative of the int rinsic time. Th e simplest choice 
for the intrinsic time flow is •& = \\dev (jValanisl . Il97ll ). However, more 
complex definitions can be given, such as: 



r) 



with ( = \\dev (i) 



(49) 



where ( is the intrinsic time scale and the positive function /i(C) = l/q(C), 

such that /i(0) = 1, is sometimes called hardening-softening function (IBazant and Bathl . 



19761). 



15 



4-1 A new formulation of endochronic models 



In this section, the endochronic model defined by Eqs. (HHj) is innovatively de- 
scribed by its Helmholtz free energy and a suitable pseudo-potential associated 
with generalized normality conditions. This approach allows for insightful com- 
parisons between endochronic models and Prandtl-Reuss models. The main 
implications will be discussed later. Let v = (e, e p , () and q nd = (cr nd , T nd , R nd ^j 
be the assumed state variables and the associated non-dissipative thermody- 
namic forces, respectively. They are the same as in the Prandtl-Reuss model 
with isotropic hardening. The Helmholtz free energy \l/ reads: 

^ = - (e - e p ) : C : (e - e p ) (50) 
2 



This form is a particular case of the one originally proposed by I Valanid (Il97ll ). 
since only one tensorial internal variable, the plastic strain, is considered here. 
The first two non-dissipative forces <r nd and T nd are the same as in Eq. (T24|) . 
while R nd = since $ is assumed to be independent of the scalar variable (. 
The pseudo-potential is defined as follows: 



|deu[C:(e — e p ) 



(i',i p ',C) e V sucli thai 
tr (i p ') =0 



(51) 







9(0 



The first term of 0, in which the stress <r nd = C : (e — e p ) is written as a 
function of the state variables, is equal to the intrinsic dissipation $ m when 
(' assumes the actual value (. The first condition associated with the closed 
convex set D introduces the plastic incompressibility assumption, while the 
second condition characterizes the plastic strain flow of endochronic theory, 
as it can be seen by comparing it to Eqs. fHHl) and fH9|) . Finally, the positivity 
of £' is imposed in order to guarantee that is positive. Using the language 
of the endochronic theory, the internal variable ( corresponds to the intrinsic 
time scale, while the intrinsic time $ is defined by its flow •& = C/g {C) ■ The 
variable ( does not directly appear in the Helmholtz free energy density and 
its associated thermodynamic forces, dissipative and non-dissipative, are thus 
zero. However, ( is not zero during the plastic evolution and plays an important 
role in the definition of i p . 
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The conjugated pseudo-potential is, in this case, of the following form: 



0* (<T d ',T d ',R d ';e,eP,() = 

= ^,,p',c)en K : ^ + ■ & + Rd 'C ~ ( 52 ) 
= I (>')+ Ie (r d ',R d ';£,eP,C) 

where E = {{r d ' , R d ') E § 2 xE such that / (V, R d '; e, e p , C) < o} and 

f (-a Rd >. e £P A _ *°(t*)Mo<* ~ £P )) \\ dev [c:(e - e^]f «,(53) 



The expression fl53|) defines the loading function of endochronic models. It is 
associated with a set E in the {r d \R d '^j space. In Figure H] this set is repre- 
sented in the case of tension-compression with g (() = 1, together with the 
projection of D on the (e p , (''J-plane. This last set is indicated by D. Some 
important remarks have to be made. First, as the system evolves, both sets 
change, due to their dependence on the internal variables. At every instan- 
taneous configurations, the set © is a straight line starting from the origin. 
The corresponding sets E are half-planes orthogonal to D. Moreover, Eq. fl50l) 
entails that R nd = —R d = and, accounting for the indicator function Io(cr d ) 
in ( |52|) . it also leads to 

T d = _ T nd = a= C . ( £ _ £ p) ( 54 ) 



Therefore, at the actual stress state (r d , R d ) the loading function / is always 
equal to zero. In other words, (r d , R d ) always belongs to <9E, during both load- 
ing and unloading phases, and all the states are plastic states. The normality 
conditions lead to the endochronic flow rules: 

<* = "t^r^ . C = A with A > (55) 

Eqs. (1521)- (|53|) and fl55|) prove that endochronic models are associative in gen- 
eralized sense. Moreover, since / is always equal to zero at the actual state, 
the loading-unloading conditions reduce to the requirement of the plastic mul- 
tiplier A to be non- negative (see the inequality in ([515]) ). In addition, the time 
derivative / at (r d , R d ), computed accounting for the fact that / also depends 
on e, e p and (, is also equal to zero and therefore, the consistency condition 
is automatically fulfilled and cannot be used to compute A. 

This situation is typical of endochronic theory and entails that the plastic 
multiplier A = ( has to be defined by an additional assumption. When the 
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function g (() is also fixed, the plastic flow i p and the intrinsic time flow 
i? = (/g (() are then known. T he standard choices are q ( () = 1 and $ = ( = 
\\dev - It has been shown in lErlicher and Pointl (120041 ) that more complex 
definitions can be chosen, such as g (£) = 1 and 



( = dev (r d ^j (l + ^sign {dev (r d ^j : e)) \dev (r d ^j : i\ 



-P < 1 < P , n>0 



(56) 



which is effectiv ely the Karray-Bouc-Casciati model (IKarray and Boud . I1989T ) 
(jCasciatil . Il989l ). It must be noticed that both flows i p and ( can be different 
from zero during unloading phases, i.e. when dev (r d ^j : i < 0. This situation, 
which is not possible in classical plasticity, occurs when 77^/?. Figure [5k 
illustrates for the mono-dimensional case the effect of n for given values of the 
other parameters: in the limit of increasing n-values the Prandtl-Reuss model 
is retrieved. Figure [5b shows unloading branches for different 7//? ratios, the 
other parameters being fixed: plastic strains may occur and tend to zero when 
7/7? tends to 1. 



4-2 Endochronic theory vs. Prandtl-Reuss model 



Consider the endochronic model, as formulated in the previous section, and the 
modified Prandtl-Reuss model. The significant state variables e, e p and £ are 
the same in both cases. Moreover, Eqs. f[3T|) and f[50l) show that the Helmholtz 
free energies differ only by the term £ ((), which is zero in endochronic theory. 
The main differences concern pseudo-potentials, as seen comparing Eqs. f[52l 
and (l5Tj) . However, the strict relationship between the two models can be 
highlighted by imposing that (' = i p in (15T1) : when (' > 0, the condition 

\\dev (C : (e — e p ))\\ = ^-g (£) must be fulfilled, while for £' = there is no 
limitation on dev (C : (e — s p )). As a result, the endochronic pseudo-potential 
([5Tp becomes equal to 



(s',s p ',C';C) = f g(C) C' + is (e',i p \C' 



i',i p ,C eV such that tr (i p =0 and C' 



(57) 



The set D and the function are not convex (see Figure [6^). However, the 
Legendre-Fenchel conjugate of is still well-posed (Appendix, item 5) and can 
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be explicitly derived from the standard procedure: 



= sup ( ^ c , )e6 (a d ' :i' + r d ' +R d '(' 
= I (V')+I E (r d ',R d ';( 



(5f 



with E = {(V, G § 2 xR such that / (V, () < o} and 



2G 



f(r d ',R d ';C) = \\dev(T d ')\\-^- g (C) + R d ' 



(59) 



Provided that ^ = J\<J y , Eqs. (1551) - (1591 also define the Legendre-Fenchel 
conjugate of the proper convex lower semi-continuous function (Appendix, 
item 5) 



4> = cl ( conv (j) 



£',£ p ',(')gV such that tr (>') = and (' > 



(60) 



which corresponds to the pseudo-potential of a modified Prandtl-Reuss model, 
in the case £ (C) = ^ = (see Eqs. (j32)). 

A similar comparison between the classical Prandtl-Reuss model (Eqs. (1231) 
and (12 5p ) and endochronic models is possible as well, but only when the former 
is perfectly plastic, i.e. if £ (£) = 0, and conditions £ (C) = and g = 1 hold 
in the latter. Note that these assumptions have been adopted in Figure [5j 



4-3 Multi-layer models of endochronic type 



The concept of assembling in parallel several plastic elements can be applied 
to the case in which each element is of endochronic type. The approach is 
analogous to the one adopted in Section 13.41 Let e and (ef , Q) be the rele- 
vant state variables. Then, the Helmholtz energy is defined as the sum of iV 
contributions, of the same kind as in Eq. (|50|) : 



N 



A' 



* = £*< = E 



i 



(61) 



i=l 



where the internal variables ef have the meaning of plastic strain of the i — th 
endochronic element. The thermodynamic forces associated with Q are zero, 
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viz. R™ d = 0. Moreover, N independent pseudo-potentials are assumed to be 
of the type (|5T|): 



|deu[C:(£ — £?) 

2G Si (Ci)/ft 



e',ef ,Ci) e V such that 



tr(4)=0, C'>0 and ef = ^ffejff] g 



(62) 



with /3j > 0, #j (Ci) > and <?j (0) = 1. The pseudo-potential of the multi-layer 
model is = J2iL± & an d its dual is (ft* = Y^iLi 4>h with 



I (cr d ')+I Ei (rf,ii!f; 



(63) 



where E { = { (rf, Rf) E § 2 xR such that /< (rf , i?f ; e, e? , 0) < o} and 

^(rf) : ^[c:(s-sf)] _ \\ devM s - sm\ 2 , Rd < (64) 

2G Si (Ci)/ft 



2G Si (Ci)/ft 



The flow rules then become of the form fl55|) . Moreover, it can be easily proved 
that at the actual state represented by (rf,Rf), the identities fi — fi — 
hold and, for this reason, the fluxes Q = Aj > cannot be computed from the 
consistency conditions and have to be defined using a further assumption. 

If the number of elements is iV = 2 , qi = q? = 1 an d both fluxes £i and C2 are 



of the form (J5b 
condition ^ = 
type. 



, then the model of ICasciatil (119891 ) is retrieved. Moreover, the 



into ( l62i) leads to a multi-layer model of Prandtl-Reuss 



5 Non-linear kinematic hardening models 



The N LK hardening rule was first suggested by lArmstrong and Frederick 
(119661 ). who introduced a dynamic recovery term in the classical Prager's linear 
kinematic hardening rule. Several modifications of this basic rule have been 
proposed, in order to improve the description of the cyclic behavior of metals. 



partic ularl y for the ratchett i ng ph enomenon (see, among others, IChaboche 



(Il99lh and lOhno and Wang! (11993h ). 
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According to traditional formulation, N LK hardening models do not fulfil 



1990l 



pp. 



the assu mption of general i zed n ormality (ILemaitre and Chaboche 
219-221) (IChaboche et al.Ul995T) . Following an approach based on the notion 
of bipotential, De Saxce ( 19921 ) introduced implicit standard materials and 



showed that the plasticity models with NLK hardening rules are of such type. 



In this section, another formulation is suggested, which leads to the proof that 
NLK hardening models belong to the class of generalized standard materials, 
provided that a suitable, non-conventional, loading function is defined. First, 
the state variables v = (e,£ p ,(, /3,(i) have to be introduced. The first three 
are the same as for Prandtl-Reuss and endochronic models, while (3 and (j 
are related to NLK hardening rule. The role of the scalar variable £i will 
be discussed later on. The corresponding thermodynamic forces are q nd = 

T d ,R d ,X. d ,Rf) e V*. The Helmholtz 



(<r nd , r nd , R nd , X nd , Rf\ and q d = (t 
energy density is chosen as follows: 



(e - e p ) : C : 



(3) : D : (e? - (3) 



(65) 



The quantity ol = e p — (3 is usually adopted as the internal variable associated 
with the kinematic hardening. However, the choice of (3 as a representative 
internal variable appears more suited, because it highlights the formal analogy 
between the first quadratic term in Eq. (|65|) . typical of plasticity models, and 
the second one, associated with the kinematic hardening. The isotropy assump- 
tion leads to the usual expression for C and entails that D =Dil ® 1+D 2 1. 
The non-dissipative forces can then be readily evaluated: 



-nd 



.nd 



X 



nd 



C : (e - e p ) 

-C : (e - e p ) + D : (e p - (3) , 
-D : (e p - (3) , 



R nd 
R nd 








(66) 



The three tensorial non-dissipative forces are related by the identity 



-nd 
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—a — X . Moreover, let the pseudo-potential be equal to 
+I 6 (i , ,& / ,(',p',&e,er,£,/3,<; 1 ) 



(i',i p ' X',(3',([) G V such that 
tr(>')=0, C'>||^'||, 

tr(/3')=0, ^=^^6, 
v ' -f si (ft) 

C( = /i(e J e p ,C,ACi)C / >0 



(67) 



with 5, g (£) , Pi (Ci) > and g (0) = (?i (0) = 1. Figure [7] shows two projections 
of the effective domain D for the tension-compression case. The first term in 
the definition of </> is identical to that of Eq. (1521 for a modified Prandtl-Reuss 
model with = 0. The second term is related to the NLK hardening and 
it is formally identical to the one used in the definition of endochronic models 
(see Eq. (I51j) ). with the substitutions dev (e) — > £ p , £ p — > (3 and C ~* Ci- The 
same analogy applies to the conditions defining the set D. 



The dual pseudo-potential then becomes 



sup/ 



e',ef ,C',/3 ,C( el 



V 



X d ' : /j' + i?f C( 



(68) 



I (<r d ') + I E (r d \ X d ', i2f ; e, e\ (, (3, Ci 



where E = { (V', X d ', i?f ) G § 2 xR x § 2 xR such that / < o} and 



f=\\dev (r d ')\\ ^ ^a y g (() + R d ' 



+ 



X d, :[D:(£P-/3) ] || D: (g-ff) ||' 



£>2 ffl(Cl)/«5 



£>2 flx(fl)/« 



+ i?f U(e,eP,C,/3,Ci: 



(69) 



Eq. (1691) defines the loading function of a model with NLK hardening and the 
associated set E is depicted in Figure [8] for the tension-compression case when 
g(Q = 1. The normality condition associated with (ft* leads to the following 
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flow rules: 



dev[T"j 
\dev( T d ) 



A = nA 

with A > 0, / < 0, Xf = 



C = A 

Ci = /i(e,e p ,C,/3,Ci)A 



(70) 



The thermodynamic force X d = — X nd = D : (s p — (3) is traceless, due to 
the assumptions adopted for the traces of i p and j3. Special attention must 
be paid to the relationship between the fluxes (\ and (. The time derivative 
of (i is defined as the product between ( and the function h, which depends 
on the state variables and must be non-negative and finite, but is otherwise 
free. The variable Ci can be interpreted as an intrinsic time scale for the NLK 
hardening flow rule. 

Accounting for the identities (r nd , R nd , X nd , Rf) = - (V, R d , X d , Rf) and 
Eqs. ( 166|) . one can prove that R d = and that the term proportional to h in 
Eq. fl69|) is always zero at the actual state. Hence, only the first two terms in 
the expression of / affect the consistency condition / = 0, which leads to the 
plastic multiplier 



(n:e) 



l + D2 _ 



1 n:X" 

2G gi(Ci)/8 



h(e,eP,CP,(i) +\{% 



2£y_ dg(Q 
3 2G dC 



(71) 



The positive functions g,gi and h determine the actual model. 

T he choice q = qi = h = 1 corres ponds to the basic NLK hardening model 



of lArmstrong and Frederick! (1l966l ). Another interesting case is given by g 
gi = 1 and 



h = 
h = 



d 2 /s 



<ki : n) if D:(£ p -/3)^0 
if D : ( £ p - (3) = 



(72) 



where mi > and ki = tt— 4 avrr is the unit vector having the same direc- 

tion as X d = D : (e p — (3). These conditions lead to 







X° 



X c 



D 2 /6* D 2 /5\D 2 /5 



(ki : e p ) = i p - 



X a 
1% 



(73) 



which is the NLK hardening rule proposed by lOhno and WangJ (1l993l ) for 
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modelling the ratchetting phenomenon in metal plasticity. It is interesting to 
compare the quantity 



<ki : i p ) 



(74) 



and the intrinsic time flow defined in Eq. f l56|) for endochronic models of the 
Bouc-Wen type. Two significant differences can be observed: (i) the governing 
flow variable is the plastic strain for NLK hardening rule and the total strain 
for the flow rule of the endochronic model; (ii) due to presence of the absolute 
value instead of the McCauley brackets, the endochronic model of Bouc-Wen 
type introduces non-zero flows during unloading phases when 7^,8. 



5. 1 From an endochronic model to a NLK hardening model 



Valanis ( 19801 ) and Watanabe and Atluri ( 1986 ) proved that a NLK harden- 
ing model can be derived from the endochronic theory by adopting a special 
intrinsic-time definition, namely when the intrinsic time scale flow ( is forced 
to be equal to the norm of the plastic strain flow. The approach suggested in 
this paper not only confirms this result, but allows for a generalization, due 
to the presence of a second intrinsic time scale (1, in general distinct from 
(. Consider the differential equations defining an endochronic model with a 
kinematic hardening variable 



-d . 



tr (&) = 3K tr (e) 

dev (&) = 2G dev(i) — j3 dev (er — 

tr (r d ) = 

T d = D 2 i p -5 T d ^-, 

z 91 (CO 

Ci = /l(£,£ p ,C/3,Cl) C 



c 

9(0 



(75) 



The idea of a kinema t ic ha rdening variable in an endochronic model was first 
suggested by iBazantl (119781 ). who however considered a linear evolution of T d 
as function of the plastic strain. An alternative way to describe the model 
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defined by (ITS]) is 



(T = C:(e-e p ) 

C = (if - §C) 1 ® 1 + 2GI 

fr(#») = 0, e^ ^-^ C 



T d = D : (e p — /3) 
D =D a l ® 1 + D 2 I 

4 ■ ( 76 ) 

Ci = A(s,£',C,ACi)C 



Moreover, both Eqs. (ITS]) and (176j) can be derived from (1651) and the following 
pseudo-potential: 



9(0 a ' ^ fli(Ci) 

+%(£',£ p, ,C , ,/3 / ,C(;^^,C,ACi 



(e'^C'^Ci) e V such that 

tr (0) =oJ= ^ Ci, Ci = A (£, e", C, A Ci) C' > o 



(77) 



Let C 



be the chosen intrinsic time definition and assume 



3 a y 

Then, introducing these conditions in (ITTI) . one obtains a pseudo-potential <fi 
which differs from the one of Eq. (IBTj) only in the inequality > s p , which 

is an equality in 0. This difference affects neither the expression of the dual 
pseudo-potential 0* = <fi* (Appendix, item 6) nor the flow rules, which are 
in both cases equal to Eqs. (I68|) - (l69|) and Eq. (ITU]) , respectively. More over, in 
the pa rticu lar case h = 1 and q ( C ) = O ] (£), the results discussed by IValanis 
( 1l980l ) and IWatanabe and Atluril (119861 ) are retrieved. 



6 Generalized plasticity models 



Generalized plasticity models (ILubliner et al.l . 119931 ) are considered an effec- 



tive alternative to NLK hardening models, since they behave s imilarly and 
are computationally less expensive (lAuricchio and Taylorl . Il995l ). A new de- 
scription of these models is suggested here, supported by a suitable pseudo- 
potential and the generalized normality assumption. In order to expose the 
basic principles of this new approach, only the simple generalized plasticity 
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model presented by lAuricchio and Taylorl (119951 ) is considered. The extension 
to more complex cases is straightforward. 



First, the state variables v = (e, e p , () have to be introduced. The correspond- 
ing thermodynamic forces are q nd = (^er nd , r nd , R nd ^j and q d = (cr d , r d , R d ^j . 
The Helmholtz energy density is chosen as follows: 



# = - (e - £ p ) : C : (e - e p ) + -e p : D : e p 

2 v ; v 1 2 



(78) 



The expression for C and D are the same as in NLK hardening models. The 
non-dissipative forces can be readily evaluated: 



r nd 



C : (e - £ p ) 



.nd 



-C : (e - e p ) + D : e p , R nd = (79) 



Note that er nd and T nd are related by the identity T nd = — (cr nd — D : e 
where the backstress D : e p introduces a linear kinematic hardening effect 
Moreover, let the pseudo-potential be equal to 



(i',i p ',('; 



e,e p 



g(e,e p ,C) (' + h (i' , i p ' , (' 



e', i p eV such that tr (i p ) = and C' > 



(80) 



where 



f(e,e*,Q ■■- 



2 

3 a y 



Hi so 

\dev[C : (e - e p ) -D : e p ] 
\dev [C : (e - e p ) -D : e*>]|| 



if /<0 
if />0 



2 

3 a y 



Hi S o C 



with H iso > 0. The main characteristic of this pseudo-potential function is 
given by the piecewise expression introduced to define the positive function g. 
It is assumed that g depends on the sign of the function /, which in turn is 
related to the state variables. The conjugated pseudo-potential 0* reads 

0* (<T d ',T d ',R d ' ] e,e p ,() 

= BUP^^ K : i' + r d ' : i p ' + R d '(' - <j>) (82) 
= Io (<T d ')+I E (r d ',R d ';e,e p ,() 
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where E = [(r d ' , R d ') G § 2 xE such that / (r d \ R d '; e, £ p , () < o} and 



f(r d \R d '; 



dev [t 
dev {r d ' 



(sjla y + H iso () + R d ' if / < (83) 

\\dev [C : (e - e p ) - D : s p }\\ + R d ' if / > 



The loading function / also has a twofold definition: recalling that the actual 
thermodynamic force T d fulfils the following identities 



C : (e - e p ) - D : e p = a - D : e p 



and R d = -R nd = 0, one can prove that if / (e, e p , () < then / (r d , R d ; e, e p , ( 
f (e, e p , C); moreover, if / (e, e p , () > 0, then / {r d , R d ; e, e p , is always zero, 
viz. the actual state represented by (r d , R d ) remains in contact with the 
loading surface dE. In Figure El this situation is depicted for the tension- 
compression case. 



The normality conditions associated with the loading function / read: 

dev(T d ^ ■ ■ 



with Xf = 



A = nA, 



C = A 



/ < A > 



These flow rules are identical to those of a Prandtl-Reuss model (see Eqs. 
(1351) ). However, they derive from a different loading function and for this 
reason the computation of the plastic multiplier A is not the same. When 
/ (r d , R d ; e, e p , (j = f (e, e p , Q < 0, the loading-unloading conditions reduce 

to A = 0, leading to an elastic behavior. As a result, the function / is also 
called yielding function, while the surface defined by the condition / = is 
called yielding surface. Conversely, when / > the set E evolves by virtue of 
the dependence of / on the state variables e, e p and (. During this evolution, 
the actual thermodynamic forces (r d , R d ) always satisfy the condition / = 0. 
Moreover, the consistency condition 



/ 



df . -d' , df hd' 
'dT dT ' + dR^ K 

_l_ M. ■ i+M. ■ iP+W C 



(86) 



,T d ,R d '=R d ) 



is also identically fulfilled and, like for the endochronic theory, it does not 
permit to compute A > 0. Hence, the condition that the so-called limit function 
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is eq ual to zero has to be invoked and this leads to (jAuricchio and Taylor 
19951): 



x = C 



1+ 



(n:S) 

N(M-f) + (D 2 +H iso ) M 



if /<0 
if < f < M 



2G f 



where M, N > 0. It can be proved that when / tends to M, the expression of 
the plastic multiplier of a classical plasticity model with linear kinematic and 
isotropic hardening is retrieved. Moreover, if Hi SO = an asymptotic value of 
||r d || exists, and is equal to \f\cFy + M. 



7 Conclusions 



A common theoretical framework between Prandtl-Reuss models and endo- 
chronic theory as well as NLK hardening and generalized plasticity models was 
constructed. All models were defined assuming generalized normality. It was 
therefore proved that a unique mathematical structure, based on the notions 
of pseudo-potential and generalized normality, was able to contain plasticity 
models traditionally formulated by other approaches. In particular, no exten- 
sion of the generalized standard class of materials had to be introduced to 
describe NLK hardening and generalized plasticity models. This approach al- 
lowed several comparisons, that have clarified the relationships and analogies 
between these, a priori different, plasticity theories. 



28 



A Appendix 



The vector spaces considered in this paper are: (i) the space of second order 
tensors; (ii) the space of symmetric second order tensors S 2 ; (iii) the set of real 
scalars R = (—00 + 00); (iv) the cartesian product of a finite number of such 
spaces. They are all equipped with an Euclidian product, so they are always 
isomorph to the Euclidian vector space X = R n . 

(1) A subset C of X is said to be: 

(a) a convex set if (1 — A) x+Ay GC whenever x, y GC and < A < 1. 

(b) a cone if Ay G C when y G C and A > 0. 

(2) Let : X — ► (—00, 00] be an extended-real- valued function defined on the 
vector space X. Then, 

(a) the epigraph of is the set 

epi — {(y, ji) such that y G X, ji G R, /i > 0(y)} (A.l) 

(b) is said to be convex on X if epi is convex as a subset of X x R. 

(c) a convex function is said to be proper if and only if the set 

© = {yGX:0(y)<+oo} (A.2) 

is not empty. The set D is called effective domain of <p, it is convex 
since is convex and is the set where is finite. 

(d) is said to be continuous relative to a set D if the restriction of 
to D is a continuous function. 

(e) is lower semicontinuous at x GX if 

(x) = lim inf (y) (A.3) 

It can be proved that the condition of lower semi-continuity of is 
equivalent to have that the lev el set {y : (y) < a} is closed in X for 



every a G R (jRockafellarl . Il969l . pg. 51). As a result, when is a proper 



convex function with a (convex) effective domain D closed in X and is 



continuous relative to D, then is lower-semicontinuous (IRockafellar 



119691 . pg. 52). 

(3) Let X* be the dual of X. Since X = R", then X** = X and the duality 
product between x and x* elements of the dual vector spaces X and X* 
can be written as x*-x. 

Let : X — ► (— 00, 00] be an extended-real-valued convex function. Then, 
the subgradients of at x G X are elements x* G X* such that 

VyGX, 0(y)-0(x)>x*-(y-x) (A.4) 
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The subdifferential set <90(x) is the set of all subgradients x* at x: 

<90(x) = {x* G X* such that the condition f[Q|) holds} (A.5) 

The function is said to be subdifferentiable at x when <90(x) is non- 
empty. 

(4) If a function : X — > (— oo, oo] is convex, proper, non-negative and such 
that 0(0) = 0, then the normality condition 

x* G <90(x) (A.6) 

viz. x* belongs to the subdifferential set of at x, entails that x*x > 0. 
Proof: Setting y = in the inequality (jA.40 entails that, for any x in 
the effective domain of 0, — 0(x) > x* • (0 — x). Hence, by virtue of the 
non- negativity of 0, x*-x > 0. 

(5) When a function : X — ► (—oo, oo] is proper, convex and lower semi- 
continuous, the dual function 0* : X* — > (— oo, oo] , defined by the 
Legendre-Fenchel transform 

Vy*GX* 0*(y*)=sup(y*-y-0(y)) (A.7) 

yex 

is related to by a one-to-one correspondence, in the sense that for such 
a kind of functions, the conju gate 0* is in turn p roper, convex and lower 
semi-continuous and 0** = (jRockafellarl . Il969l . pg. 104). 
Under these assumptions, it also holds: 

Vy*GX* 0*(y*)=sup(y*-y-0(y)) (A.8) 

yGD 

Moreover, the following relationships are equivalent: 

(i) x*G90(x); 

(ii) x G <90* (x*) 

(iii) 0(x) + 0*(x*) = x*-x (A.9) 

Condition (i) is equivalent to x* • x— 0(x) > x* • y— 0(y). The supremum 
of the second term of this inequality is equal by definition to 0*(x*) and 
occurs when y = x and therefore (iii) is the same as (i). Dually, (ii) and 
(iii) are equivalent. 

Remark 1. Under the previous assumptions, if > and 0(0) = 0, then 
f)A.7j) entails that 0*(O) = 0. Moreover, the identity 0** = implies that 
0(0) = sup y * eX , (— 0*(y*)), which in turn leads to 0* > 0. Reciprocally, 
0* > and 0*(O) = entail that > and 0(0) = 0. 
Remark 2. If 0* is such that 0* > and 0*(O) = 0, then the normality 
condition (ii) implies that x* • x > 0. Proof: Condition (ii) is equivalent 
to (i), with > and 0(0) = 0. Then, using the result of item 4, the 
non-negativity of x* • x follows. 

Remark 3. The conjugate 0* of an arbitrary function : X — > (— oo, oo] 
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can still be defined by (1 A. 71) . In this case, 0* is proper, convex, lower semi- 
continuous and is equal to the conjugated 0* of = cl (conv <pj , where 
is the greatest proper c onvex lower semi-continuous function majorized 
by (f> ( iRockaferlari . Il969l . pp. 52, 103-104). 
(6) A function : X — > (— oo, oo] is positively homogeneous of order 1 if and 
only if 



VyGX, VpG (0,Oo), 0(py)=p0(y) 



(A.10) 



The epigraph of such functions is a cone (iRockafellarl . Il969l . pg. 30 ) . 
Given : X — > (— oo, oo], the following three statements are equivalent: 

(i) is proper, convex, lower semi-continuous and positively homoge- 
neous of order 1 

(ii) The Legendre-Fenchel conjugate 0* of is the indicator function of 
a non-empty, convex and closed set E, i.e. 



</>* (y*)=W) 



if y* G E 
-oo if y* G E 



(iii) is the support function of a non-empty, convex and closed set E, 
i.e. 

0( y ) = I|(y) = sup (y*-y) 



The equivalence between (i) and (i i) can be proved b y showing that 0* 
has no values other than and +oo ( IRockafellarl . Il969l . pg. 114). The set 
where 0* = is non-empty, convex and closed since is proper, convex 
and lower semi-continuous. The equivalence between (ii) and (iii) follows 
from the definition of Legendre-Fenchel transform, support functions and 
indicator functions. 

Remark. If fulfils conditions in (i), then for any x where is subd- 
ifferentiable, 



(x) = x*-x with x* G <90(x) 



Proof: from the equivalence between (i) and (ii), the conjugated of is 
the indicator function of a closed convex set E and x* G E since is 
subdifferentiable at x by assumption. Then, use Eq. flA.9[) and recall by 
(ii) that 0* (x*) = 0. 
(7) Let : X — ► (— oo, +oo] be a proper, convex, lower semi-continuous 
function, positively homogeneous of order 1. Then: 
(i) from item (6), its conjugate 0* is the indicator function of a non- 
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empty, closed and convex set E. Hence, by using the definition (1A.4I) . 



dh fx*) 



if x* G int ( E 

C (x*) if x* G dE 
if x* i E 



(A.11) 



where C (x*) = jx G X : Vy* G E x- (y* - x*) < o} is the so-called 
normal cone at x* G <9E. 
(ii) if in addition does not depend on some components yi of y = 
(yi,ya) C X = Xi x X 2 , i.e. 0(y) = 0(yi,y 2 ) = <j>(yz), then the 
conjugated function 0* can be computed as follows: 



'(y*. y*2) = su P( yi ,y 2 )ex(y^ ■ yi + yl ■ y2 - 0(y 2 )) 

= Io(yJ) + sup y2eX2 (y^ ■ y 2 - 0(y 2 )) 

= I (y 1 )+I E (yl) 



(A.12) 



The Legendre-Fenchel conjugate is the indicator function of with 
respect to y* plus the Legendre-Fenchel conjugate of 0(x 2 ), which is 
the indicator function of a non-empty, closed and convex set E. Hence, 

x G dig (x*) ^> xi G Xi and x 2 G 9I E (x 2 ) 

In the particular case where E= {y* G X* such that / (y*) < 0} , 
where / is a convex and smooth function, the normality condition at 
y* = x*, viz. x G 8Ie (x*), can be written as follows 



x = \i grad f (x*) 



with 



H = for / (x*) < 
H > for / (x*) = 



These two last conditions are often replaced by 
/i > 0, / (x*) < 0, fif (x*) = 



(A.13) 



which are the classical loading-unloading conditions of plasticity, usu- 
ally written with /x replaced by the plastic multiplier X. The dependence 
of / on the argument x* is often omitted in order to simplify the nota- 
tion. In the convex mathematical program ming literat u re. (|A.13D are 
known as Kuhn- Tucker conditions (see e.g. iLuenbergerl (119841 ) ) . 
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Fig. 1. Classical Prandtl-Reuss model. Tension-compression case. 

a) Projection of the pseudo-potential effective domain D on the (i p , ('^J -plane. This 

set is indicated by D. b) Domain E associated with the dual pseudo-potential <f>* . 
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Fig. 2. Modified Prandtl-Reuss model. Tension-compression case with = 0. 

a) Projection of the pseudo-potential effective domain D on the (^i p , £M -plane. This 

set is indicated by B. b) Different configurations of the domain E. The position of 
E changes according to the value of the internal variable £. The point (T d ,R d ), 
representing the actual state, always lies on the axis R d ' = 0. 
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Fig. 3. Mechanical dissipation for the case of simple tension. The hatched area is 
the energy J <f> m (t)dt dissipated during the monotonic loading, 
a) Classical Prandtl-Reuss model, b) Modified Prandtl-Reuss model, c) Modified 
Prandtl-Reuss model with £(£) = ^|(£) = 0. d) Modified Prandtl-Reuss model with 

V§3c(0 = a y(9(0 ~ 1) : t ne classical model is recovered. 
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Fig. 4. Endochronic model. Tension-compression case with g(Q = 1. 

a) Several configurations of the set B, which is the projection of the pseudo-potential 

effective domain D on the (i p -plane, b) Configurations of the convex set E 

associated with those of ID>. The point (r d , R d ), representing the actual state, always 
lies on the axis R d = 0. 
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Fig. 5. Endo chronic Karray-Bouc-Casciati model (thin lines) vs. Prandtl-Reuss 
model (thick line). Tension-compression case with g(Q=l and £(£) = 0. 
a) Influence of the parameter n on loading branches, b) Influence of the 7//? ratio 
on unloading branches. The slope at an = is the same for all 7//? values. 
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Fig. 6. Endochronic model vs. Prandtl-Reuss model. Tension-compression case, 
a) The set ID is the projection of on the (i p , £'J -plane, where B is the non-convex 

effective domain of the pseudo-potential <p of Eq. (|57p . It defines an endochronic 
model where the intrinsic time flow £ equals the norm of i p . b) The convex set 
E associated with the indicator function (j)* given in Eqs. (|58p -(|59 p . which is the 
Legendre-Fenchel conjugated of (j>. 
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Fig. 7. NLK hardening model. Tension-compression case. 

a) Projection of the effective domain D on the (i p , -plane, b) Projection of D 
on the m', C{ J-plane. 
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Fig. 9. Generalized plasticity. Tension-compression case. 

a) Projection of the pseudo-potential effective domain D on the (i p , £M -plane. This 

set is indicated by B. b) Several configurations of the domain E. When / > 0, E 
translates upward during loading phases and downward during unloading phases. 
The point (T d ,R d ), representing the actual state, always lies on the axis R d ' = 0. 
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